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We study the Lagrangian mechanism of the fluctuation dynamo at zero Prandtl number and 
infinite magnetic Reynolds number, in the Kazantsev-Kraichnan model of white-noise advection. 
With a rough velocity field corresponding to a turbulent inert ial-range, flux-freezing holds only in a 
stochastic sense. We show that field-lines arriving to the same point which were initially separated 
by many resistive lengths are important to the dynamo. Magnetic vectors of the seed field that 
point parallel to the initial separation vector arrive anti-correlated and produce an "anti-dynamo" 
effect. We also study the problem of "magnetic induction" of a spatially uniform seed field. We 
find no essential distinction between this process and fluctuation dynamo, both producing the same 
growth-rates and small-scale magnetic correlations. In the regime of very rough velocity fields where 
fluctuation dynamo fails, we obtain the induced magnetic energy spectra. We use these results to 
evaluate theories proposed for magnetic spectra in laboratory experiments of turbulent induction. 
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I. INTRODUCTION 

Hannes Alfven introduced the notion of magnetic flux 
conservation in ideal magnetohydrodynamics [1 . This 
property of "flux-freezing" involves an implicit assump- 
tion, however, that plasma fluid velocities remain smooth 
in the limit of vanishing resistivity. This assumption need 
not be valid if the viscosity vanishes together with the 
resistivity, at a fixed or a decreasing magnetic Prandtl 
number. At the implied high Reynolds numbers, smooth 
laminar plasma flows will be unstable to the development 
of turbulence. We have argued previously that magnetic 
flux-conservation in its usual sense cannot hold at any 
small resistivities, no matter how tiny, in a turbulent 
plasma with a Kolmogorov-like energy spectrum [2] . Nei- 
ther is flux-freezing completely broken but, instead, it 
becomes intrinsically stochastic due to the roughness of 
the advecting velocity field [3] H] . Infinitely-many mag- 
netic field lines are carried to each point — even in the 
limit of very high conductivity — which must be averaged 
to obtain the resultant magnetic field at that point. 

The previous ideas are theoretical predictions for tur- 
bulent flows of astrophysical and laboratory plasmas. 
However, they are exact properties of at least one sol- 
uble problem of turbulent advection, the Kazantsev- 
Kraichnan (KK) model of kinematic magnetic dynamo. 
See [5HZ] and many following papers [5HT6]. In that 
model the physical turbulent velocity field is replaced 
by a Gaussian random field that is white-noise in time 
and rough (Holder continuous) in space. It is rigorously 
established for this model that Lagrangian particle tra- 
jectories are "spontaneously stochastic" in the limit of 
infinite Reynolds numbers, with an infinite ensemble of 
trajectories for a given advecting velocity field and a 
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given initial particle position [T7H24] . This remarkable- 
breakdown of Laplacian determinism is a consequence of 
the properties of two-particle Richardson diffusion [25] , 
which are also expected to hold for real fluid turbulence. 

The stochastic nature of flux-freezing must play an es- 
sential role in the turbulent magnetic dynamo at high ki- 
netic Reynolds numbers, either for fixed or for vanishing 
magnetic Prandtl numbers. The zero Prandtl- number 
fluctuation dynamo was studied in an earlier work of 
Eyink and Neto [26], hereafter referred to as I, primar- 
ily for the KK model. This model undergoes a remark- 
able transition as the scaling exponent £ of the velocity 
structure function is decreased, with dynamo effect dis- 
appearing for £ < 1 in space dimension three. This seems 
rather counterintuitive, because velocity- gradients and 
line-stretching rates in the model increase as £ is lowered, 
which should favor dynamo action. Paper I showed that 
the stochasticity of flux-freezing is fundamental to under- 
stand this transition. Although magnetic line-stretching 
rates indeed increase for decreasing £, the resultant mag- 
netic field at a point is the average of infinitely- many field 
lines advected to that point from a large spatial volume. 
As shown in I, the correlations between independent line- 
vectors are small for £ < 1, when the velocity field is too 
rough, and the net magnetic field decays despite the rapid 
growth in strength of individual field lines. 

The contribution of flux-line stochasticity to turbulent 
dynamo action for £ > 1 was not fully explicated in I, 
however. Infinitely-many field lines are advected to each 
point from a large spatial volume, e.g. with diameter 
of order the velocity integral length L in one turbulent 
turnover time L/u rms . But all of these field lines are un- 
likely to make an equal contribution to dynamo action. 
In particular, field vectors that arrive to the same point 
from very distant initial locations must be poorly cor- 
related and give small (or negative) dynamo effect. It 
is not hard to guess that dynamo action in the limit of 
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high- Reynolds number Lu rms /u ^> 1 and low Prandtl 
number ^ « 1 must arise mainly from vectors with 
initial separation distances r of order ~ the resistive 
length. Indeed, this is the only available length-scale in 
that limit, on dimensional grounds. The exact contri- 
bution to fluctuation dynamo from vectors with initial 
separation r remains unknown, however, even in the KK 
model. One of the main purposes of this paper to provide 
a quantitative answer to this problem. 



Another important issue left unresolved in I was the 
outcome of a spatially-uniform initial magnetic field in 
the dynamo regime of the KK model for £ > 1. Can 
such a magnetic field provide the seed for a small-scale 
fluctuation-dynamo? This question is very closely con- 
nected with the proper formulation of the concept of 
"magnetic induction" . This process is usually defined 
as the creation of small-scale magnetic fluctuations by 
the turbulent "shredding" of a non-zero mean magnetic 
field and has been invoked to explain small-scale mag- 
netic fluctuations in liquid metal experiments [27-31], in 
related numerical simulations [32, 33 , and at the solar 
surface [34j [35] . It is frequently suggested that this pro- 
cess is distinct from fluctuation dynamo. For example, 
Schekochihin et al. wrote: "Given a multiscale observed 
or simulated magnetic field, one does not generally have 
enough information (or understanding), to tell whether 
it has originated from the fluctuation dynamo, from the 
mean-field dynamo plus the turbulent induction or from 
some combination of the two." [33] However, if a mean 
magnetic field can provide a seed for fluctuation dynamo, 
is there really any meaningful physical distinction be- 
tween the two mechanisms? Or if magnetic induction 
is defined instead as a process that occurs only in the 
absence of fluctuation dynamo, as it sometimes is, then 
how can one regard the two processes as acting in "com- 
bination" ? Clearly there is some conceptual confusion in 
the literature about the proper definition of magnetic in- 
duction. We shall use our exact results in the KK model 
to discuss these issues and, also, to evaluate some the 
proposed theories of the induced magnetic spectrum. 



The contents of this paper are as follows. In section [TT| 
we provide more detailed background to our work and 
a more formal statement of the problems to be stud- 
ied. In section [TTT] we study the Lagrangian mechanism of 
the fluctuation dynamo at zero Prandtl number and the 
quantitative role of stochastic flux-freezing. A first sub- 
section |III A| presents exact mathematical analysis and 
the second subsection IIIB gives numerical results. The 
next section [IV] discusses the problem of magnetic induc- 
tion using our analytical results for the KK model. We 
draw our final conclusions in section [V| Appendices [A] 
and [B] present some technical material. 



II. BACKGROUND AND PROBLEM 
STATEMENT 

We begin this section with a summary of the key results 
of Paper I. A "dynamo order-parameter" was introduced 
there with purely geometric significance. Let x(a, t) be a 
stochastic Lagrangian trajectory, satisfying 

dx = u(x,f)cft+ y/2rjdW(t), x(t ) = a. (1) 

We shall assume in this section that the advecting veloc- 
ity field is incompressible, V*u = 0, as in I. W(t) in 
is a vector Brownian motion and 77 is the resistivity (or, 
in fact, the magnetic diffusivity ) . Let ^(a, t) be a pas- 
sive vector starting as a unit vector at space point a 
and time to and subsequently transported along x(a, £), 
stretched and rotated by the velocity- gradient. That is, 



4(a,t) = e fc -V x(a,t). 
Then I introduced the quantity 



(2) 



K ke (r,t) = (4(MK(a',^ 3 (x(a,t)-x'(a^))), (3) 

with r = a — a 7 , where (•) denotes average over the tur- 
bulent velocity realizations, (•) denotes average over the 
random Brownian motions and t indicates a second, inde- 
pendent realization of the latter. This quantity measures 
not only the magnitude but also the angular correla- 
tion of independent line- vectors that start as unit vectors 
£fc(a,to) = e /c? -^( a '^o) = &t at points a, a 7 separated by 
displacement r, which end at the same point at time t. 

Now let B be a magnetic field kinematically trans- 
ported by the random velocity field u and diffused by 
resistivity 77: 

ftB + (u.V)B = (B.V)u + 77AB, B(t ) = B (0) . (4) 

As shown in I, the mean energy in this magnetic field can 
be expressed as 



{B\t)) = J c 



d 3 r^(r,t)(4 0) (r)5f 0) (0)). (5) 
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It has been assumed in deriving ^ that initial condi- 
tions on the magnetic field are statistically independent 
of the velocity and that both quantities are statistically 
homogeneous (space translation-invariant). This formula 
makes clear the close relation of dynamo action to the 
line-correlation (3j), which grows exponentially rapidly 
precisely in the dynamo regime. 

This connection was verified in I by an exact calcula- 
tion of the line-correlation (|3| in the soluble Kazantsev- 
Kraichnan (KK) model of turbulent dynamo [5H7J fTUI 
[T6] . This model replaces the true turbulent velocity by a 
Gaussian random field which is white-noise in time: 



(^(x,tV (x',0) = i^ ij (r)5(t - t f ) 



(6) 



with r = x — x'. In this model, the magnetic correlation 
function 



C ij (r,t) = (B i (r,t)B^(0,t)) 



(7) 
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evolves according to a closed, linear equation 

d t C ij = M%C U , (8) 

where M is a singular diffusion operator. If K % i (r) ~ 
for small r, with rugosity exponent < £ < 2, then the 
velocity realizations are rough in space and provide a 
qualitatively good model of turbulent advection. It was 
shown in I that 7£/^(r, t) in the KK model is the solution 
of the adjoint equation 



The quantity (12) measures the correlation of two mde- 



(9) 



with initial condition 7£^(r,£) = 5keS 3 (r). This fact was 
exploited to show that 7£^(r, t) in the KK model at zero 
magnetic Prandtl number, as expected, grows exponen- 
tially in the dynamo regime (£ > 1) but has only power- 
law time-dependence in the non-dynamo regime (£ < 1). 

The study in I left open, however, several important 
questions. First, what range of r = |r| substantially con- 
tributes to the correlation 7Z and to the space-integral 
in (J5|? It is known in the KK model that magnetic 
field lines arrive to a point from a spatial region with 
radius L(t) ~ (t — to) 1 ^ 2- ^, m the root-mean-square 
sense. This is also expected for kinematic dynamo in a 
real turbulent fluid, where £ = 4/3 and L(t) ~ (t-t ) 3/2 
corresponds to turbulent Richardson diffusion of the field 
lines. However, not all of the field lines in this large vol- 
ume should be expected to contribute equally to dynamo 
action. Field vectors starting from points separated by 
distances r much larger than the resistive length-scale 
may arrive at the same point with small angular cor- 
relations and thus contribute little to net magnetic field 
growth. It follows from the analysis of I for the KK model 
in its dynamo regime that, for times t long compared to 
the resistive time-scale, 



Hki(r,t) ~e- xt g M (r), (A < 0) 



(10) 



where —A is the largest (positive) eigenvalue of M* and 
Q is the corresponding (right) eigenmode. Thus, ques- 
tions about the spatial structure of 1Z reduce, in the ex- 
ponential growth regime at sufficiently long times, to the 
corresponding questions about the right eigenmode Q of 
M* (or, equivalently, left eigenmode of M). It is one 
of the goals of the present work to calculate this eigen- 
mode and thereby determine the degree of correlation of 
line- vectors arriving from various initial separations. 

A second important issue left unresolved in I was the 
response of a spatially uniform initial magnetic field to 
turbulent advection and stretching. In the case of a uni- 
form initial field, formula (|5| simplifies to 



(B 2 (t)}=n ke (t){B^Blo)), 



(11) 



with 



n M (t) = J d 3 rK M (v,t) = (4(£)<W)o. ( 12 ) 



pendent line- vectors that arrive at the same point, re- 
gardless of their initial separation. The result (10) does 



not imply exponential growth of this integrated correla- 
tion, however, unless one can show that 



d 3 rg k£ (r)^0. 



(13) 



Below we shall demonstrate this fact. 

This result is closely connected with another impor- 
tant physical problem, the creation of small-scale mag- 
netic fluctuations by turbulent "magnetic induction" of 
a non-zero mean-field [27H33] . Indeed, the formula (11) 
applies directly to this situation, which corresponds to 
the special case (B^B^) = (B^)(B^). This observa- 
tion already makes clear that there is no essential physi- 
cal distinction between "fluctuation dynamo" and "mag- 
netic induction", which simply correspond to different 
choices of initial magnetic field (random vs. determin- 
istic). Other definitions of "magnetic induction" are of- 
fered in the literature. If the fluctuation dynamo fails to 
operate for some reason (e.g. if the magnetic Reynolds 
number is too small) , then the term "magnetic induction" 
is sometimes used for the process by which small-scale 
fluctuations are generated from the mean magnetic field. 
If a mean-field dynamo still operates and creates an ex- 
ponentially growing large-scale magnetic field, then these 
"parasitic" small-scale fluctuations may be also exponen- 
tially increasing. This provides another mechanism to 
produce exponential growth of small-scale magnetic fields 
in the absence of a fluctuation dynamo. If instead there 
is no mean-field dynamo, then such growth of small-scale 
fluctuations by "magnetic induction" is sub-exponential. 
We shall discuss below the non-dynamo regime (£ < 1) 
of the KK model as a simple example of such "magnetic 
induction" and use it to evaluate some physical theories 
that have been proposed for the phenomenon. 

We mention finally one additional motivation to study 
the left eigenmodes of M.. As pointed out in I, these have 
a physical interpretation as correlators of the magnetic 
vector potential A. More precisely, the correlation func- 
tion defined by 



g ke (r,t) = (A fc (r,t)Aj(0,f)), 
in the KK model satisfies 



d t Q, 



kt 



(14) 



(15) 



with diffusion operator adjoint to that in eq. ((§. The un- 
derlying reason for this fact is the conservation of mag- 
netic helicity H = f d 3 r A«B for the ideal induction 
equation. If the latter is written as 



d t B = £B, 



(16) 



for a linear operator £ (depending on the velocity u), 
then helicity conservation is equivalent to 



d t A = -C*A, 



(17) 
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up to addition of a gradient term. Different operators C 
are possible, which coincide in their action on solenoidal 
magnetic fields (V«B = 0), but whose adjoints £* cor- 
respond to different gauge choices for A. Equation (17) 
directly implies fll5|) i n the KK model. The addition of 
resistivity to eqs. Q16|), (17) does not change this result, 
since the additional Laplacian operator is self-adjoint. 



III. FLUCTUATION DYNAMO 

We study in this section the Lagrangian mechanism of 
the fluctuation dynamo in the KK model for three space- 
dimensions (3D) and for homogeneous and isotropic 
statistics of both velocity and magnetic fields. We be- 
gin in the first subsection |III A with exact mathematical 
analysis of the problem. Most of this discussion applies 
to a fairly general situation, allowing for compressibility 
and helicity of the velocity field. We then specialize to the 
case of an incompressible, non-helical advecting velocity, 
with a power-law space-correlation corresponding to an 
infinite- Reynolds-number inertial-range range for the ve- 



locity field. In the second subsection IIIB we present nu- 
merical results for this latter case and discuss their phys- 
ical interpretation. The reader who is mostly interested 
in the final results, and not their detailed derivation, may 
skip directly to that section of the paper. 



A. Mathematical Analysis 

We discuss two natural choices for the linear operator 
C in the ideal induction eq.(16) and the corresponding 
adjoint operators C* and gauge choices for the vector 



potential, successively in the following two subsections. 



1. Gauge I 



One choice to write the induction equation is as 



d t B = Vx(uxB -77VXB). 



(18) 
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Note that in the ideal equation with 77 
— Vx(uxB) is the Lie-derivative acting on B as a dif- 
ferential 2-form. See [36 , eq.(2.3) for VB = 0. Thus, 
this form of the induction equation is most directly con- 
nected with the conservation of magnetic flux through 
advected 2-dimensional surfaces. The corresponding ad- 
joint equation for the vector potential is 



d t A = ux(VxA)-r/Vx(VxA) 
= -(u-V)A + (VA)u 

+r/[AA- V(V-A)]. 



(19) 



In the gauge-choice implied by this equation, pure-gauge 
fields A = VA satisfy V(c^A) = and are, thus, time- 
independent up to possible spatially constant terms. 



Boldyrev, Cattaneo and Rosner in [14] developed an el- 
egant formulation of the KK model based on the eq.(18) 
for the case of homogeneous and isotropic statistics. We 
review their formulation — hereafter referred to as the 
BCR formalism — in Appendix [Aj paying special atten- 
tion to issues of gauge-invariance. Readers not previously 
familiar with the work of [14] may wish to review that 
appendix before proceeding. One essential observation of 
BCR is that the evolution operator for magnetic correla- 
tions factorizes as M = VJ. More precisely, in the KK 
model starting from eq.(18) 



where 



d t C ij = M%C pq = D ijM JkL Pq C pq 



V ij,k£ = e ikp e j£qQ Q 



(20) 



(21) 



is the non-positive, self- adjoint differential operator 
which relates the magnetic correlation to the vector- 
potential correlation and 



ji 



ij,k£ ^ikptjiq 



(22) 



is the self-adjoint multiplication operator with T pq (r) = 
2r]S pq + K pq (0) - K pq (r), where K™(r) is the spatial 
velocity-correlation in eq.Q. Explicitly, 

M%C pq = d r d s (T rs C ij ) - d p d s (T is C pj ) 

-d r d q (T rj C iq ) + d p d q (T ij C pq ) . (23) 

One may further write equations in the KK model for 
the joint correlations of magnetic and vector-potentials 

¥ fc (r,t) = {B\v,t)A k {0,t)), * fc *(r,t) = {A k {v,t)B\0,t)) 



using the operator 1Z % 



ipk dr>, in terms of which the 



operator V itself factorizes as V^ ,ki 
d t W n = -{K i ' m )*J mn , Pq 



n \* q 1? p > k \b q 

<J rnn/pq /v * 



= lZ^ k W^. Then 
K q ^\ 

k 



(24) 



Finally, the vector-potential correlation in ([14]) obeys 

d t Gk£ = (M*)%G P 
with M* = JV. Explicitly, 



Jk£,rs^ TS,Pq Qp 



(25) 



(M*)%g pq = T rs d r d s g M - T ps d k d s g pi 



-T rq d r d £ g kq + T pq d k d £ g pq . (26) 

BCR mainly considered the isotropic sector of the KK 
model, invariant under proper rotations but not neces- 
sarily space reflections. In that case, with r = r/r, 

K lj (r) = K L (r)r l r j + K N (r)(5 lj - r l r j ) + ^H(r)e ljk r k , 

where kl, kn, and kh — rg are conventional longitu- 
dinal, transverse and helical correlation functions, with 
similar representations of C, (5, etc. However, ref.[l4] in- 
troduced instead a special basis to expand tensor corre- 
lation functions, as C ZJ = J2a=i CaCdf w ^ n 



V2r 



(5 l - 



fV), # 



S3 
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so that 

d = y/2rC N , C 2 =rC Ll C 3 = V2rC H . 

The special feature of this basis is that the Hilbert-space 
inner-product defined by 

(G,C) = J d 3 r ^(r)C^(r) 

simplifies in the isotropic sector to 

(G, C) = 4tt / dr [G1C1 + G 2 C 2 + G 3 C 3 ]. 
./o 

We hereafter use the notation C, (5 etc. in the BCR for- 
malism for the isotropic sector to represent the column 
vectors 




In this representation, the operators J and V take the 
simple forms 




with 

a(r) = V2[2rj + k n (0) - K N (r)] 
b(r) = 2r] + « L (0) - « L (r) 

c(r) = V^(0)-#(r)]r. 

A crucial observation of BCR is that V factorizes as 
V -711V with 

/0 d r °\ / 

ft = ^ \,K*=\-d r & 
\0 ~^d r r 2 J V r 2 d r ^ 

Another important fact is that Ran(7£), the range of the 
operator 7£, consists of the solenoidal functions C ZJ that 
satisfy c^C* J = djC 1 ^ = 0. Thus, all solenoidal solutions 
of the equation dfC — M.C = —7Z7Z* JC can be obtained 
by solving 

d t W = -1VJ1IW (27) 
with a self- adjoint operator <S = and then setting 

d r W 2 

-^d r (r 2 W 3 ) 



C = KW = 



Comparison of eq.(27) in the isotropic sector with the 
general eq.([24|) suggests that W is a simple linear trans- 
formation of the magnetic-field, vector-potential correla- 
tion For this result, and for the derivation of all the 
preceding statements, see Appendix |A) 

We now discuss the solutions of the adjoint problem 
dtG = M*G = — JVSW Q . Its solutions are likewise re- 



lated to solutions of (27) by the equation 





W = 1l*G 



-d r G 1 + ^G 2 
r 2 d r 



This relation is many-to-one. Since Ker(7£*) = 
[Ran (72.)] , the kernel of 7Z* consists of the correlations 
of gradient type: 



Gi(r) = V^A'(r), G 2 (r) = rA"(r), 



(28) 



for a scalar correlation function A. Thus, any solutions 
Q of the adjoint equation that differ by such a gradient 
solution are mapped to the same solution W of (27). This 



freedom just corresponds to gauge- invar iance of W. 

We next employ this formalism to discuss the right 
eigenfunctions C a and left eigenfunctions Ga of which 
satisfy 



MC a 



"Aq-Cq,, M*Got — —^aGa, 



respectively. (Of course, for the continuous spectrum of 
the operators these are generalized eigenfunctions.) Let 
us suppose that one has obtained the eigenfunction W a 
of the self- adjoint operator S = TV J1Z: 



sw a = n*jnw a = \ a w a . 

Then it is not hard to see that 

1 



C a =KW a , Ga 



A, 



-jnw a 



(29) 



(30) 



The first result is obvious. To obtain the second, note 
that, with Ga as defined above, 

n*Ga = ^-sw a = w a . 

^a 

If this result is used to eliminate W a in the definition 
of Ga 1 then the statement follows. An interesting con- 
sequence is that right and left eigenfunctions are related 
very simply by 



Ga 



1 

An, 



\J &a • 



Note also that {Cp,Q a ) (7lWp,j^JKW a ) 

^(WfrSWa) = (WfrWa) = 5 ai p. Thus" the left and 
right eigenfunctions are biorthogonal. 

Our discussion so far has been fairly general, per- 
mitting a compressible velocity field with reflection-non- 
symmetric (helical) statistics. However, we now special- 
ize. If space-reflection symmetric statistics are assumed, 
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then W 3 = C 3 = G 3 = 0. Furthermore, c(r) = in the 
operator J. Thus, the eigenvalue problem ( |29| ) reduces 
to a single equation for W 2 : 



-d r {b{r)d r W 2 ) 



a(r) - ra'(r)]W 2 = XW 2 . 



This is a standard Sturm-Liouville eigenvalue problem. 
The relations ( 30 ) yield 



Ci = d r W 2 , C 2 



V2 



w 2 



(31) 



(true even in the reflection non-symmetric case) and 



Gi 



V2a(r) 



Xr 



W 2 



b(r) 



d r W 2 , G 2 



a(r) 



d r W 2 . (32) 



The problem further simplifies if one assumes an incom- 
pressible flow, implying the relation a = \f2 (b + \rb') . 
In that case the Sturm-Liouville problem becomes 



- d r (p(r)d r W 2 ) + q(r)W 2 = XW 2 . 

with p(r) = b{r) and 

_ 2b(r) - 2rb'(r) - r 2 b"(r) 
Qv) — o — —u r 



(33) 



;d r (r 2 b(r)) 



From general Sturm-Liouville theory, the spectrum is 
bounded below and may consist of both point spectrum 
and continuous spectrum, depending upon the choice of 
the function b(r). Dynamo effect corresponds to the low- 
est eigenvalue being negative, A < 0, with a ground state, 
square-integrable wavefunction drW^ir) < oo. The 
ground-state eigenfunction W 2 > if it exists at all, must be 
positive for all r > by the Sturm oscillation theorem. 



2. Gauge II 

There is another natural choice for the linear operator 
C in the ideal induction eq.(16). This choice corresponds 
to writing that equation as 



d t B = -(u-V)B + (B-V)u - (V-u)B + rjAB. (34) 



If B is a smooth solution of eq. ([34]) for 77 = and 
if p is the mass density solving the continuity equation, 
dtp + V-(pu) = 0, then, as is well-known, i = ~B/p is a 
"frozen-in" field. That is, 



_9 

dt 



i = -(u-V)£+ (^V)u. 



The righthand side of this equation is just C^£, the 
Lie-derivative operator acting upon vectors (rank-1 con- 
travariant tensors) [36]. These facts explain the interest 
of the eq. (34), since intuitive geometric notions of ma- 



terial line-vector motion can be exploited to understand 



magnetic dynamo effect. The results in I all depended 
upon using this form of the induction equation. 



The equation for the vector-potential adjoint to (34) is 



d t A = -(u- V)A - (Vu) A + 77AA. 



(35) 



We use the tilde to distinguish the vector potential in 
this gauge from that resulting from (19). There is 

also a geometric significance to eq.(35), since £u^A = 
(u-V)A + (Vu)A is the Lie-derivative acting on A as 
a differential 1-form. This is related to the fact that 
"frozen-in" magnetic flux through surfaces (2-cells) can 
be written as the line-integral of the vector potential A 
(or A) around closed loops (1-cycles). In the gauge choice 
corresponding to (35) the pure-gauge fields A = VA have 
zero Lagrangian time-derivative, D t A = 0, up to possible 
spatial constants. Of course, the two gauge choices A and 
A must be related by a suitable gauge transformation 

A = A - VA. 



In fact, eq.(35) can be rewritten as 



with 



d t A = ux(VxA) - t?Vx(VxA) - V$, 



$ = u-A — 7?V-A, 



implying that A = / dt' $(£')• 

In the KK model, the form of the induction equation 



( 34 ) leads to the evolution equation for the magnetic cor- 



relation d t C lJ = M l p J q C pq , with 



M%C pq = d r d s (T rs C ij ) - d s (d p T is C pj ) 

-d r (d q T rj C iq ) + (d p d q T ij )C pq . (36) 



Comparison with the definition of M. in ([23]) shows that 
MC — MC when acting on elements C of the subspace 
of solenoidal correlations functions. We note further- 
more that this subspace is invariant under the evolu- 
tion (J8| (because the induction equation (18) preserves 
solenoidality of B.) The adjoint equation for the vector- 
potential correlation Qij(r,t) = (Ai(r,i)Aj(Q,i)) is 



d t G, 



Wk£ 



(37) 



with 



T rs d r d s g k£ + (d k T ps )d s g p£ ^ 
Hd £ T rq )d r g kq + (d k d £ T pq )g pq . (38) 



In the previous section we have discussed how to deter- 
mine the eigenvalues and eigenfunctions of the operator 
M.. The right eigenfunctions C a of M. and C a of M. are 
the same in the solenoidal subspace, since M. = M. there. 
We shall obtain the left eigenfunctions of M. by solving 
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for the differences with the corresponding left eigenfunc- 
tions ofM. Defining bQ a =Q a -Qa, JM* =M*-M*, 
it follows from M*Q a = -^ a G a ,M*G a = -X a G a that 



( y M*+\ a ) sg a = -SM*g a . 



(39) 



A straightforward calculation using ([26]) and (38) gives with 



(5M*)l q e G Pq = d k [T pq (d q g pi - d £ g pq )] 

+d £ [T pq (d p g kq - d k g pq )] 
+d k d £ (r pq g pq ) . (40) 



The solvability condition of eq.(|39| is (C a ,SM*g a ) = 0, 
for the right eigenfunction C a of M (which is also the 
right eigenfunction C a of M) with eigenvalue X a . This 
is easily seen to be satisfied using the definition of 5M.. 
Thus, eq.(39) has a unique solution 5g a in the subspace 
orthogonal to C a . Defining Q a = g a + Sg a , it follows that 



Similar arguments using eq.(|39|) show that 5M*g a and 
5 g (x are orthogonal to every solenoidal eigenfunction 
C(3 — Cp and, thus, to the entire subspace of solenoidal 
correlation functions. Hence the new set of eigenfunc- 
tions C a , g a for all a form another biorthogonal set. 

There are several simplifications in the isotropic sec- 
tor of the model. As shown in Appendix |Aj any func- 
tion Sg in the isotropic sector which is orthogonal to all 
solenoidal correlations must be pure-gauge: 



Sg k £ = d k diA. 



(41) 



(For simplicity, we drop here the spectral index a which 
labels the eigenvalues and eigenfunctions.) Finding 
Sg thus reduces to finding A. Note furthermore that 
A4*5g = when Sg is pure-gauge (ultimately because 
gauge functions A do not evolve in the gauge-choice of 
eq.([T9|)). Thus, the equation for Sg becomes 

(5M* + \) 5g = -5M*g. 



and, substituting from (41), 



d k d e (T pq d p d q A) + xd k d £ A = -dk [T pq (d q g pi - d £ g pq )] 

-d £ [T pq (d p g kq - d k g pq )\ 
-d k d £ (T pq g pq ) . 

Finally, there must exist in the isotropic sector a scalar 
function $(r) such that 



t™ (d p g kq - d k g pq ) = a fc $, 



(42) 



with d k <&(r) = fk^(r), so that the previous equation 
becomes 

T pq d p d q A + AA = - [2$ + T pq g pq ] . 
Substituting the isotropic form 

T pq (r) = T L (r)r p r q +T N (r) (5 pq - r p r q ) +T H (r)6 pqm r m , 



and the similar expression for g into (42) gives, after 



some computation, $'(r) = 2(T N ^ H -T H ^ N ). The re- 
sulting equation to be solved for A is 



T L d 2 r A + 2T N -d r A + AA = - [2$ + T pq g pq ] , (43) 



$( r ) 



dp [T N (p)* H (p) - T h (p)*n(p)] 



(so that $(+oo) = 0) and 

T pq G pq = T L G L + 2T N G N + 2T H G H . 



Alternatively, the equation (43) may be written using 



BCR quantities, with Tl — T/v = a/y/2 
# = aW2 ~ cW * - V2c* L , 



*i(r) = V2 [ dp 
Jo 



\ W 3 (p) 



and 



TG = T pq G n 



ad + bG 2 + cG 3 



B. Numerical Results and Physical Discussion 

As a particular case of physical interest we shall con- 
sider an incompressible, non- helical velocity field with 
a = y/2(b + \rb') and c = 0. To model an infinite- 
Reynolds-number inertial-range, we take, for < £ < 2, 

b(r) = 2r) + 2D^. 

We non-dimensionalize the equations of the KK model 
using £rj = (rj/Di) 1 ^ for space and = for time. 
The Sturm-Liouville problem is then given by (33), or 



-d r {p{r)d r W 2 ) +q(r)W 2 = XW 2 , 



p(r) =2(l + r*), 



4 2(£ + 2)(£-l) 



r 2-£ 



(44) 



Endpoints r = 0, oc are singular, non-oscillatory and 
limit-point. The Friedrichs boundary conditions to se- 
lect the principal solution are W^(0) = ^2(00) = 0. E.g. 
see [37]. The Sturm-Liouville operator with q(r) above 
is obviously non-negative for £ < 1, so dynamo effect 
can exist only for £ > 1. In the latter case, there is both 
point spectrum and, above some threshold value, A > A c , 
continuous spectrum. The lowest eigenvalue is negative, 
implying dynamo effect. To give the closest correspon- 
dence with real hydrodynamic turbulence, we take the 
Richardson value £ = 4/3 in our work here. See [12] for 
a numerical study of general values 1 < £ < 2. 
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FIG. 1. Plot of ground-state eigenfunction W2(r) versus r. 
All quantities in this and following figures have been non- 
dimensionalized with resistive units as discussed in the text. 



FIG. 2. Small-r asymptotics of ground- state eigenfunction. 
In blue is a log-log plot of W2<V) versus r and in red the 
quadratic fit (0.1288)r 2 . 



We obtain the solution of (|44j) using the MATSLISE soft 
ware for MATLAB 38] [39 



Briefly, this package solves regu- 
lar Sturm-Liouville eigenvalue problems by first convert- 
ing them through the so-called Liouville transformation 
W2 = ^/[p(r)] 1 / 4 and x = J Q r -^==into a correspond- 
ing Schrodinger equation for ^(x). The latter is then 
solved by a high-order Constant Perturbation Method 
(CPM) algorithm, which generates its own numerical grid 
of points r^, i = 1, 2, L. To treat our singular problem 
on an infinite interval, we truncate to a finite interval [0, b] 
with the condition W2(b) = and then take successively 
larger b values, as discussed in [39] , Section 6.3.1. We 
have considered successively b = 200, 700, 800, with 
convergence of solutions for r < 500. The ground- state 
eigenvalue found by this method is A = —0.1927, in good 
agreement with [12], Fig. 5, and the corresponding eigen- 
function W 2 (r) is plotted in Fig. 1. The eigenfunction is 
normalized so that 



f 

Jo 




drWi(r) = 1. 



(45) 



FIG. 3. Large-r asymptotics of ground- state eigenfunction. 
In blue is a log-linear plot of ^(r) versus r 1 ^ 3 and in red a 
straight line with slope from (47) for £ = 4/3. 



Asymptotics at small and large r are known for the 
ground-state eigenfunction. The small-r behavior was 
obtained by [13], who noted 



W 2 (r) At 2 + Br 2 ^ + o(r 2+ ^), r«l. 



(46) 



Indeed, direct substitution of this ansatz shows that the 



Sturm-Liouville eigenvalue equation (44) is then satisfied 
up to terms o(r^), if and only if B = —A. The large-r 
behavior from a WKB asymptotic analysis is p~2j [13] 



W2M ~ exp 



V^2Ar( 2 ^)/ 2 

2-e 



r > 1 



(47) 



up to a power-law prefactor. 
the dominant balance 2r^W2 ' 



This corresponds to 
— XW2 in the Sturm- 
Liouville eq.(|44j) at large-r. Both the small-r and large- 



r asymptotic behaviors predicted analytically have been 
verified in our numerical solution. As shown in Fig. 2, the 
leading-order r 2 behavior in (46) is verified over about 24 



orders of magnitude with the value A = 0.1288 obtained 
by taking the small-r limit of the numerical results for 
W2/r 2 . Although we shall not show it here, we have fur- 
thermore verified the subleading r 10 / 3 term in (46) with 



the same value of A, over about 16 orders of magnitude. 
We also verify the stretched-exponential decay (47) at 
large-r, as shown by the log-linear plot of W2(r) vs. 



,1/3 
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FIG. 4. Plot of magnetic field correlation functions, longitu- 
dinal Cl{t) in blue and transverse Cn{t) in versus r. 




-G L 



FIG. 5. Plot of magnetic vector-potential correlation func- 
tions, longitudinal Gl(t) in blue and transverse Gn(t) in 
green versus r. 



in Fig. 3. The red line shows the slope predicted by (47) 
with A = —0.1927. Our numerical evaluation of W 2 {t)i^ 
in very good agreement with all known analytical results 
for the dynamo eigenfunction. 

Magnetic and vector-potential correlations in the dy- 
namo growth mode are obtained from W 2 via (31), (32) 
and plotted in Figs. 4 and 5 using the traditional lon- 
gitudinal and transverse functions. The normalizations 
are those which follow from (45). One obvious feature 



is the much longer range of the vector-potential correla- 
tions as compared with the magnetic-field correlations. 



The small-r behavior is obtained from (|46j) to be 

C L ~ y/2A[l - + o(r% (48) 



C N ~ V2A[1 - (1 + |)r* + o(r% (49) 



and 



G L ~G 



N 



4\/2A 
A 



[1 + 0(7 



(50) 



Note that the contributions of order cancel in the latter 
two functions, signalling their greater smoothness. The 
large-r behavior is found, with (47), to be 



A /2 
r-Z /2 W 2 , C 2 ~ —W 2 
r 



2 + t* 



(51) 



(52) 



The signs implied for r ^> 1 are Cjv < 0, Cx > and 
Gat, Gl > 0. Although it is difficult to see clearly in 
Fig. 4, Cn{t) < in our numerical solution for r > 6.85. 
The existence of negative tails in Cat was noted some time 
ago [H [9] to be a consequence of the solenoidal character 
of the magnetic field (see also just below). 

The physical-space behaviors discussed above imply 
corresponding results for the magnetic energy spectrum 
of the dynamo mode. At low wavenumbers 



E(k) ~ A'k A , k<^k v 



(53) 



with krj = 2ir/t rt and A' some constant numerically pro- 
portional to A. To see this, note that the stretched- 
exponential decay (51) implies that E(k) is analytic at 



low-/c and can be expanded as a convergent power- series 
in k 2 . The leading term proportional to k 2 vanishes, how- 
ever, since the integral J d 3 rC^(r) = 0. By isotropy, this 
is equivalent to the vanishing of the integral J dr 
^ Cl+C2 is the trace C T 



r 2 C T , 



where Ct 



Cu. Using (31), 



2 C T = V2(rd r W 2 + W 2 ) = V2d r (rW 2 ) , 



so that 



POO 

/ drr 2 C T (r) = V2 lim rW 2 (r) = 0. 
Jo ' r ^°° 

Note that this result requires the existence of negative 
tails in Ct at large r [8, 9 . The k A spectrum in (53) was 
explained long ago by Kraichnan and Nagarajan[6] as 
due to the dipole magnetic field that results from "irreg- 
ularly twisted and elongated current loops whose trans- 
verse dimension is ~ /c" 1 [ our fc" 1 ]." 

The spectrum at high wavenumbers is [13] 



E(k) ~Cfe _(1+0 , k^krj. 



(54) 



10 




< 

I 10- 4 - 

^ 10- 6 - 



-numerics 
-(11.93)r 2 



50 100 150 200 250 300 350 400 



FIG. 6. Plot of gauge-change function A(r) versus r. 



This follows mathemati cally by Fourier transforming the 
singular terms ex in (48), (49). The result (54) is the 



exact analogue for the Kazantsev model of the Golitsyn- 
Moffatt /c -11 / 3 spectrum [40l|4l]. The dominant balance 
in the induction equation for length-scales t <C is 

d t B - r]AB = B^Vu, 

with the magnetic field at scale £ v . The time- 
derivative must be included because of the rapid change 
in time of the velocity field. Solving for the statistical 
steady-state of the above linear Langevin equation leads 

(B 2 ) 

to E(k) rsj ±—±E u (k), where E u (k) is the energy SpeC- 
^B 2 ) 

trum of the velocity field, and not to E(k) ~ ^ r ^E u {k) 1 
as in the original argument of Golitsyn. Our discussion 
here exactly parallels that of Frisch and Wirth [42] for the 
analogous passive scalar problem. Note finally that the 
energy spectrum in the Kazantsev model, both in its low- 
k and high-/c behaviors, is thus close to that argued for 
kinematic dynamo in inertial-range hydrodynamic tur- 
bulence by Kraichnan and Nagarajan [6]. 

We shall now turn to the evaluation of the gauge-II 
correlations. For this purpose, we must determine the 



gauge-change function A(r) which solves eq.(43): 



r- 2 (r 2 6A 7 ) 7 + AA 



(55) 



with F = 2<3> + TG. We have solved this equation numer- 
ically, by the procedure described in Appendix [B] The 
result is plotted in Fig. 6. 

At small r, the solution A = A reg + A sing is the sum 
of a regular part 



A 



reg 



A + Air 2 +o(r 2 ) 



and a singular part 

A sina ~A s r 2+ « + o(r 2+ «). 



(56) 



(57) 



FIG. 7. Small-r asymptotics of gauge-change function. In 
blue is a log- log plot of A(r) — A versus r and in red the 
quadratic fit (11.93)r 2 . 



-numerics 
-asymptotic; 




FIG. 8. Large-r asymptotics of gauge-change function. In 
blue is a log-linear plot of |A(r)| versus r 1//3 and in red a 
straight line with slope from (60 ) for £ = 4/3. 



This is verified by substituting into eq.(55), using the 
similar expansion for F 



F - F + Fir^ 



which follows from 



Jo P 



with $o = - J °° dr < and 

TG= 8 ^ [ 3 +(3 + 0r* + o(r*)]. 
A 
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Thus, F = 2$ + 2 ^£ A < and F 1 = ^^(3 + f) < 0. 
Then it is straightforward to obtain 

AA + 12Ai = -F (58) 



2 [2Ai + (2 + A a ] (3 + = -Fi . (59) 

These results have been verified in our numerical solu- 
tion, with specific values of the constants 



Ao 



-1535.5, Ai = 11.9299, A s 



-6.0232. 



See Appendix [Bj and Fig. 7 for a comparison of A(r) and 
A reg {r) at small r. 



l2A r (2-0/2 



2-£ 



up 



At large-r we expect A(r) ~ exp 
to power-law prefactors. Rewriting the eq.([55j) for A as 

bA" + r~ 2 {r 2 b)'A f + AA = —F 

one finds that the terms bA" ', A A cancel to leading order. 
For r > 1. 



TG< 



2 ( 2 + r 3£/2- 



^2, 



and $' = aW 2 /r - V^(2 + O^" 1 W 2 , so that 
^>^- 2(2+ ^ ) r ^/ 2 - 1 ^ 2 . 



It follows that 



F = 2$ + TG - 



The dominant balance of r 2 {r 2 b)' A' and — F gives 



A' 



x/2 

X 1 



r > 1. 



(60) 



This asymptotics is verified in our numerical solution. 
See Fig. 8 for a log-linear plot of |A(r)| vs. r 1 / 3 , with the 
red line having the slope predicted by (60), (47). 

From the function A we obtain SQ via eq.(41). We 



plot in Fig. 9 the longitudinal and transverse components 
5Gl and 5Gn obtained from (28). It is interesting that 



these functions show almost exactly opposite behaviors 
of algebraic signs compared with Gl, Gat in Fig. 5. Both 
5Gl and SGn are proportional to at small r because 
of the A s contribution and are stretched-exponentials 



^ /2 W 2 , SG 2 



~V2 



W 2 (61) 



at large r, so that SGn > 0, 5Gl < for r ^> 1. 

Finally, we obtain Q = Q + SQ. Plotted in Fig. 10 are 
Gl, Gat. These are the central results of this paper, be- 
cause of their relation to 1Z in eq.(10). We can thus 



interpret these functions as the correlations in the dy- 
namo growth phase of vector line-elements carried to the 
same final point that were initially separated by distance 




FIG. 9. Plot of change in vector-potential correlation func- 
tions, longitudinal 5Gl(t) in blue and transverse SGN(r) in 
green, versus r. 




FIG. 10. Plot of line- vector correlation functions, longitudinal 
Gl(t) in blue and transverse Gn(t) in green, versus r. 



r, with directions longitudinal and transverse to the sep- 
aration vector, respectively. We shall discuss their physi- 
cal interpretation below; first we consider the asymptotic 
behaviors for small-r and large-r. Combining results for 
Q and SQ at r <C 1 gives 



G L 



' Ay/2 A 



■2Ai +(2 + 0(1 + 0A.r« (62) 



and 



G 



N 



^^+2A 1 )+(2 + 0A s r«. 



(63) 
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FIG. 11. Contribution of transverse seed field vectors to a 
positive line- vector correlation. Two transverse field- vectors 
on parallel lines at the initial time (left), when brought to the 
same point by a fluid motion, are positively correlated at first 
contact of the lines (right). 



FIG. 12. Contribution of longitudinal seed field vectors to 
a negative line- vector correlation. Two longitudinal field- 
vectors on the same line — or two adjacent lines — at the initial 
time (left), when brought to the same point by a fluid motion, 
are negatively correlated at first contact of the lines (right). 



It is interesting to note that (5, unlike (5, contains terms 
oc and is singular at small r. The reason for this is 
that the eq.(35) for A contains the velocity-gradient Vu, 



whereas the eq.(19) for A contains only u itself. Finally, 



comparing the large-r asymptotics in (52) with that in 
(61 ) gives for r ^> 1 



G 1 



G 2 



V2 



W 2 , 



(64) 



so that G N > 0, G L < for r > 1. 

The results in Fig. 10 quantify the significance of mag- 
netic field-line stochasticity for the small-Pr m turbulent 
dynamo. The plotted correlations represent the contribu- 
tion to magnetic energy at a given point produced by field 
vectors that arrive from points separated by distance r 
initially. As one can see, the contribution is quite diffuse 
(in units of the resistive length £ v ) with fat, stretched- 
exponential tails. Both correlations are positive for small 
separations r < 3^. Particularly interesting is the long 
negative tail for the longitudinal correlation Gl, which 
represents an "anti-dynamo effect" that suppresses the 
turbulent growth of magnetic field energy. 

It is easy to understand the positive signs of Gl, Gn 
at small r values. For r < £ n the Brownian motion 
term dominates in the equation which follows from (HI) 
for the relative motion r(t) = ^(t) — xi(t) of a pair 
of points. Since the Brownian motion corresponds to a 
simple translation of the lines, with no rotation, the line- 
vectors which start parallel (either longitudinal or trans- 
verse) at separation r < £ v will tend to remain nearly 
parallel when they intersect. Hence, they arrive at the 
final point positively correlated. 



The signs of the tails of Gl, Gn for large r are less 
easy to understand, because turbulent advection effects 
become important for r > £ r] . We suggest here a heuris- 
tic explanation. In the left half of Fig. 11 we show a 
pair of parallel field-lines separated by r with their ini- 
tially transverse line-vectors indicated in red. We then 
show on the right the result of a large-scale fluid motion 
which brings the two line-vectors together to the same 
point, stretched and rotated by the flow. (It is impor- 
tant to understand that the velocity at scales r > £ v 
cannot actually bring the two lines to touch, but only 
within distance ~ £ r] of each other; it is then the action 
of two independent Brownian motions which completes 
their transport to the same point.) Around the time of 
first contact shown in the figure, the two line- vectors are 
positively correlated. Further rotation by the flow can 
eventually produce ant i- alignment and negative correla- 
tion. The situation is exactly the opposite, however, for 
the initially longitudinal line- vectors shown in Fig. 12, 
which start on the same (or nearly the same) field-line. 
As these vectors are brought together by a similar fluid 
motion as before, the line vectors are anti-aligned near 
the time of first contact. They may then become posi- 
tively oriented as further rotation and stretching occurs. 
These geometric considerations may help to explain why, 
for large enough r, positive alignment dominates in Gn 
while negative alignment dominates in Gl- Transverse 
line- vectors seem to be brought together from large dis- 
tances by mainly lateral advection of field-lines, while 
longitudinal vectors are brought together by twisting and 
looping of field-lines. 
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IV. UNIFORM FIELDS AND MAGNETIC 
INDUCTION 

We now turn to the second question left open in I, 
the fate of spatially-uniform magnetic seed fields in the 
dynamo regime of the Kazantsev model. This issue is 
closely related to the problem of "magnetic induction", 
which is often defined simply as the generation of small- 
scale magnetic fluctuations by turbulent tangling of the 
lines of a non-vanishing mean-field. In the case of ho- 
mogeneous statistics — which we only consider in this 
paper — any non-zero average field must be spatially uni- 
form. Thus, the response of a uniform field to homoge- 
neous turbulence is intimately related to that of magnetic 
induction. In fact, we shall argue that there is no real 
distinction in this context between "magnetic induction" 
and "fluctuation dynamo" with a uniform seed field. 

Let us consider then the special case that the initial 
magnetic-field B( ) is spatially-uniform (but still possi- 
bly random). Applying eq.(ll) gives the energy in the 

(B)to be 



fluctuation field b = B 

(b 2 (t)) 



= Tl M {t){b k {0) b{ Q) ) + [R M (t) - d ke }(Bf 0) )(Bl 0) ). (65) 

In the second line of the equation above, one may make 
a formal distinction between "fluctuation dynamo" rep- 
resented by the first term and "magnetic induction" rep- 
resented by the second term. However, their asymptotic 
growth rates at long times are the same, both given by 



n M (t) 



(A<0) 



Qkt = I d 3 r Gki(r) 



if the integral 



is non- vanishing. We show now that it is indeed non-zero. 

Since the above integral is gauge- invariant, one may 
just as well consider the integral of Gki( r )- Furthermore, 
as long as the turbulent velocity field (but not necessarily 
the magnetic field) is statistically isotropic, then 

Qkt J d 3 r Gki(r) ^G T S ki 

where G T = 4tt / °° dr r 2 G T {r) and G T {r) = G L (r) + 
2Gn(t) i s the trace function. From ([32]) 



G T (r) 



V2d + G 2 _ 2aW 2 + (a + V2b)rd r W 2 
r Xr 2 



Using the incompressibility condition a = y/2 (b + \rb f } 
and integrating by parts once in r gives 



drr 2 G T (r) = — ^ , 

V ZA Jo 



(r 2 b" + 4rb')W 2 (r) dr. 



On the other hand, multiplying the eigenvalue equation 
( 33 ) by r 2 and integrating by parts in r twice gives 



/>00 f'OC 

- / (r 2 b" + Arb')W 2 (r) dr = X r 2 W 2 (r) dr. 
Jo Jo 



We finally obtain that 

G T = 4tt / drr 2 G T (r) = 2nV2 / r 2 W 2 (r) dr > 0. 
Jo Jo 

We used here the positivity of the ground-state (dynamo) 
eigenfunction W 2 . Notice that we have not assumed any 
specific form of 6(r), as long as the dynamo mode exists. 

In fact, the entire spatial structure of small-scale fluc- 
tuations for "magnetic induction" is the same as that 
for "fluctuation dynamo" at long times. Both are deter- 
mined by the same dynamo eigenmode. To see this, we 
note that the magnetic correlation can be expanded, in 
general, into right eigenmodes of A4 : 



(66) 



The expansion coefficients are determined from the initial 
correlation C^L (r) by 



d 3 r SajijOOCgjtr), 



(67) 



using biorthogonality of right and left eigenmodes C a and 

the 

(0)^(0)) = 



Q a . If C(Q)( r ) = (^(o)^(o)) ^ s r - m dependent, then for the 



ground state mode (say a = 0) Co = Go-ki(B^B 
lG T (Bf 0) ) > 0. Thus, with r fixed, 



c e- A °^'(r), 



t —> oo. 



The magnetic correlation C* J (r,£) is dominated by the 
leading dynamo eigenmode Cq (r) with subleading con- 
tributions from additional eigenmodes (point spectrum 
of M) if any. There is no essential physical differ- 
ence between the two cases of "fluctuation dynamo" 
(B( ) = b(o)) and "magnetic induction" (B( ) = (B( ))). 

It is interesting to consider in (66) also the opposite 
limit of large r with t fixed. In that case, the result with 
initial condition C^(r) = (B l (0) B J (Q) ) is 



'(0)^(0)/ 



^(r,t)~(BJ 0) ^ 0) 



oo. 



(68) 



This is clearly true for the case C^(r) 



<^ )X^0)). 



because the above limit then corresponds to "clustering" 
of the magnetic correlation function. However, the ex- 
pansion coefficients (67) are linear in the initial correla- 



tion, so the above result must hold whenever C^(r) is 



r-independent. The limit (68) comes entirely from the 
continuous spectrum of M in the expansion ( [66] ) (for 



which the sum over a is actually an integral), because the 
finite number of dynamo eigenmodes C a associated to the 
point spectrum of M all have the stretched-exponential 
decay in (47), with the corresponding value of X a < 0. 
Their contribution thus vanishes in the large-r limit. 

We conclude that there is no essential distinction in the 
dynamo regime between "magnetic induction" for a uni- 
form mean-field and "fluctuation dynamo" for a homo- 
geneous (but random) seed field, since all of their physi- 
cal behaviors and underlying mechanisms are exactly the 
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same. Another possible definition of "magnetic induc- 
tion" which is employed in the literature is the growth 
of small-scale fluctuation fields generated from a non- 
vanishing mean-field, but only in a parameter regime 
where the "fluctuation dynamo" does not exist. The 
term is used this way in discussion of liquid-metal lab- 
oratory experiments operated at a sub-threshold regime 
[28[ l29[ I3T] . Although such inductive growth of fluctua- 
tions is generally sub-exponential, it may be exponential 
if the average magnetic field is itself exponentially grow- 
ing because of a mean-field dynamo effect. Of course, 
with this definition, it makes no sense to regard "mag- 
netic induction" and "fluctuation dynamo" as possible 
co-existing and competing mechanisms for small-scale 
magnetic field growth. 

We may consider as an example of this second type 
of "magnetic induction" the failed-dynamo regime of the 
KK model for £ < 1, at zero Prandtl number and in- 
finite magnetic Reynolds number. Unlike the liquid- 
metal experiments where fluctuation-dynamo fails be- 

(c) 

cause Rem < Re m , the failure here is due to the extreme 
roughness of the advecting velocity field (see I). This 
model problem is not a perfect analogue of the induc- 
tion phenenoma seen in the experiments. For example, 
we solve the KK model with velocity integral scale L = 00 
(and thus Re m = 00) so that the large-scales of our prob- 
lem are quite different than the non-universal, inhomo- 
geneous and anisotropic conditions seen in experiments. 
Also, the average field is uniform in the KK model with 
homogeneous statistics. Since spatially constant fields 
must be time-independent, we cannot study in this set- 
ting induction of an exponentially growing mean-field. 
On the other hand, the cited experiments [28] [29] [31] 
have studied the turbulent induction of a near-uniform, 
stationary external field. The analogy with those exper- 
iments is close enough that we can test some proposed 
theories of "magnetic induction" in our model situation. 

We therefore consider in this light several results of I 
for the failed dynamo regime of the KK model with £ < 1. 
It was found there for the case of spatially-uniform and 
isotropic initial data C l ^(r) = A5 lJ that, at long times, 



Id 1/7 



(69) 



with 



Ci 



2 2 -§v 1 -l*« +2 > 



which is negative and decreasing from to —2 for < 
£ < 1. Thus, the energy in the magnetic fluctuations 
is growing, but only as a power-law in t. These same 
results hold , in fact, for general uniform initial data of 



the form (r) 



(BfoB 3 , Q ^). This follows directly from 



eq.(ll) when the random velocity (but not necessarily 
the magnetic field) is statistically isotropic and the line- 
correlation matrix 1Zki(t) = ^7Z(t)5k£. 

Paper I showed further that there are three spatial 
regimes of the magnetic correlation. These are the re- 



sistive range: 



1-2 



r<^, (70) 



the quasi-steady, inertial-convective range: 



C L (r,t)~A' 



L(t) 



Ci 



In < r < L(t), 



(71) 



and the very large-scale range: 



C L (r,t) ~A 



1 - 



2Ci(Ci+3 + 



my 



r > L(t). 



1 \ r 

(72) 

We have here defined L(t) = {Dit) 1 ^ 1 to be a character- 
istic large length-scale of the magnetic field. Note that 
A f ,A" are constants numerically proportional to A. Re- 



sult (70) follows from the formula for T in (a) on p. 26 of 
I, together with the series expansion of the hypergeomet- 
ric function 2 i*i (eq.(75) in I). Equations fn} and f72| 



follow from I, eq.(67) for a = with p<l and p^> 1, re- 
spectively. For the last case, the large- argument asymp- 
totics of the Kummer function in [43 , eq. 6.13.1(2), is 
employed to derive the 1/r 7 term. 

These three ranges all have simple physical descrip- 
tions, which are easiest to discuss based on the corre- 
sponding magnetic energy spectra obtained by Fourier 
transform. In the resistive range the result is 



E(k,t) oc 



(B 2 (t)} 
V 



k > k v . (73) 



(Note that r] = D\l^.) This is a Golitsyn-like spectrum, 
with physics exactly like that discussed earlier for the 
resistive range of the dynamo growth modes. This range 
is very universal in the small-Pr m limit. 

The inertial-convective range has energy spectrum 

E(k,t) ~ A'mt))^^- 1 k L (t) < k < k v (74) 

with &£,(£) = 2tv /L(t) — ^ as t — ^ 00. This growing 
range is responsible for the energy increase in (69). It is 



"quasi-steady" in the sense that all its statistical charac- 
teristics are identical to those in the forced steady-state 
(see [10] and I), with the uniform initial field replacing 
the role of the external force. The physics involves a 
nontrivial competition of stretching B*Vu and nonlinear 
cascade u-Vb. If one assumes that "induction" of the ini- 
tial, background field B( )-Vu dominates the stretching, 
then this is the same physics invoked by Ruzmaikin and 
Shukurov [44 for induced fluctuations. Assuming a bal- 
ance of the terms (u-V)b ^ b-Vu ~ B^-Vu, they pro- 
posed on dimensional grounds an Ak~ x spectrum, with 
A oc (B^ s). However, such a spectrum only occurs in 
the KK model for £ = 0. In that case, the cascade term 
u-Vb can be represented as an eddy-diffusivity 2Di and 
a balance equation 

d t h - 2L>iAb = B (0 )-Vu 
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predicts the correct spectrum (B 2 ^)k 1 . There is no dis- 
tinction in this case between the inertial-convective range 
and the resistive range with Golitsyn spectrum. For 
< £ < 1, however, there is an "anomalous dimension" 
|Ci| which corrects the dimensional prediction. There is 
no longer any simple, heuristic argument for the magnetic 
spectral exponent in the quasi-steady range, which re- 
quires the computation of a non-perturbative zero- mode. 
The physics involves a nontrivial competition between 
creation of magnetic energy by stretching and destruc- 
tion by turbulent cascade to the resistive range. 

For the very low wavenumber range k <C k^t) one 
might expect the dimensional argument of Ruzmaikin 
and Shukurov to work after all and to correctly yield an 
Ak~ x spectrum. The results in (|7Q|)-(72) correspond to a 



This leads to the prediction that 



"(0) 



)E(k) 



member of the one-parameter family of self-similar decay 
solutions found in I, with parameter choice a = 0. For 
the general member of this family, the low-wavenumber 
spectrum is of the form Ak 01-1 ("permanence of large 
eddies"; see I), so setting a = should naively lead 
to the Ruzmaikin- Shukurov prediction. But this is not 
the case. Fourier transformation of the initial uniform 
magnetic field instead yields a term in the energy spec- 
trum oc AS(k), a delta-function at k = 0. The actual 
energy spectrum in the low-wavenumber range is found 
by Fourier transforming the second term in ( 72 ) to be 



E(k,t) ~ A l [L(t)} 1 k 1 ~ 



k<^k L (t). (75) 



The simple physics of this range is direct induction from 
the spatially-constant field via the balance 



B( )-Vu. 



Indeed, solving this linear Langevin equation yields 
E(k,t) ex (B^)(Dit)k 2 E u (k), reproducing the above 
spectrum. Note that [L^t)] 1 = D\t corresponds to diffu- 
sive spectral growth in this range, due to the white- noise 
in time character of the advecting velocity field. The 
total energy in this range remains constant in time, how- 
ever, because shrinking of the range exactly compensates 
for increase of the energy spectrum. Once the energy 
b\ c± kE(k, t) in an interval around wavenumber k ap- 
proaches the energy (B 2 Q ^) in the initial uniform field a 
nonlinear cascade begins and that wavenumber k joins 
the inertial-convective range with spectrum (74). 

The above argument resembles one which has been in- 
voked to explain a k~ 5 ^ 3 spectrum of magnetic energy re- 
ported at wavenumbers k < k^ in some low-Pr m liquid- 
metal experiments [29J E] and in a related numerical 
simulation [32]. These works study the induction of an 
imposed magnetic field by a turbulent flow with a Kol- 
mogorov energy spectrum. The argument made by those 
authors amounts to assuming a balance between induc- 
tion of the external magnetic field B( ) and convection by 
the large-scale velocity U( ) (both mean and fluctuating 
components): 



(B 2 Q ^)E u (k), so that E u (k) ~ e 2/ ' 6 k 5/3 implies a similar 



r 2/3^-5/3 

magnetic energy spectrum E(k). We find this argument 
unconvincing. In the first place, large-scale advection 
conserves the energy in magnetic fluctuations. It cannot 
therefore balance the input from induction of B( ) . If non- 
linear cascade can be ignored, then magnetic energy must 
be expected to grow, as in our result ( 75 ) above, and not 



to saturate. Furthermore, if the large-scale sweeping were 
important, then it should appear in our balance argument 
for (75), because the KK model contains such dynamical 



u 



(0) 



Vb 



B( )'Vu. 



sweeping effects. However, we see that the correct spec- 
trum at very low wavenumbers is obtained by ignoring 
such sweeping as irrelevant. This is consistent with the 
results of Frisch and Wirth [42] for the diffusive range 
of a passive scalar, who emphasized the (nontrivial) fact 
that sweeping by large-scales plays no role in the balance 
for that high- wave number range. Finally, we note that 
£,-5/3 ^ e (Obukhov-Corrsin) spectrum expected for a 
passive scalar in a Kolmogorov inertial-range and not for 
a passive magnetic field. The argument of [29j [31] [32] 
thus omits the effects of the small-scale stretching inter- 
action (b'V)u and we know of no good justification for 
doing so in a steady-state, saturated regime. 



V. CONCLUSIONS 

The main purpose of this paper was to quantify the 
importance of stochasticity of flux-line freezing to the 
zero-Prandtl- number turbulent dynamo. It is a rigor- 
ous result for the Kazantsev-Kraichnan model that La- 
grangian trajectories becomes intrinsically stochastic in 
the inertial range with velocity roughness exponent £, due 
to Richardson 2-particle turbulent diffusion. In this sit- 
uation, infinitely-many magnetic lines in the initial seed 
field at time to are brought to a point at time t from a re- 
gion of size L(t) rsj (t — to) l ^ 2 ~^ . Note that the size of the 
region sampled is independent of the resistivity. Not all of 
the magnetic field lines in this large region will, however, 
make an equal contribution to the net dynamo growth 
of magnetic energy. The contribution from lines initially 
separated by distance r is quantified, in isotropic non- 
helical turbulence, by the line-correlations R^r^t) and 
-R/v At long times, these correlations are propor- 
tional to e~ xt GL(r) and e -tA£ GAr(r), respectively, where 
Gl(t) and Gn{t) are longitudinal and transverse com- 
ponents of the (left /adjoint) dynamo eigenfunction. The 
central results of this paper are these two functions, plot- 
ted in Fig. 10, which quantify the relative contribution 
to magnetic energy from lines initially separated by the 
distance r, in resistive units. 

The principal contribution to the dynamo growth 
comes from lines at separations of the order of magni- 
tude of the resistive length £ r? . However, the decay of the 
correlations, in units of the resistive length, is a slow, 
stretched exponential. This implies that lines which ar- 
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rive from any arbitrary, fixed separation r — no matter 
how large — will contribute an amount of energy growing 
exponentially rapidly in time. Of course lines separated 
by r much, much greater than £ v will make a small rel- 
ative contribution, but even lines separated by many £ v 
make a substantial contribution. To underline this fact, 
we plot in Fig. 13 below the eigenfunctions Gl{t) and 
Gn(t) multiplied by r 2 , which, integrated over r, give 
the mean magnetic energy for a uniform seed field. (In 
order to make this plot, we have extended our numer- 
ical results to r > 163^ with the asymptotic formulas 



(64).) As this figure should make clear, lines separated 



by many hundreds of resistive lengths are important to 
the dynamo growth. To be more quantitative, we find 
that lines initially separated by up to 968^ must be con- 
sidered in order to get 90% of the total magnetic energy. 
Another feature dramatically illustrated in Fig. 13 is the 
strong anti- dynamo effect arising from lines- vectors ini- 
tially parallel to the separation vector r, represented by 
the long negative tail in v 2 Gl{t). We have proposed a 
heuristic explanation for this interesting effect in terms 
of twisting and looping of field lines (Figs. 11 and 12). 




FIG. 13. Contributions of line-correlations to the energy. The 
integrands, r 2 Gl(t) in blue and r 2 Gn{t) in gre e n, fo r mag- 
netic energy with a uniform seed field in eqs.( 11 ),( 12 ). 



As we shall see in a following paper [45 these principle 
conclusions apply not only in the Kazantsev-Kraichnan 
model with Pr m = but also to kinematic dynamo in 
real hydrodynamic turbulence with Pr m ~ 1. 

A second purpose of this paper was to discuss the 
meaningful distinction, if any, between "fluctuation dy- 
namo" and "magnetic induction." In the case of the KK 
model at Pr rn = and Re rn = oo, we have shown that 
a uniform mean field (or a uniform random field) may 
provide the seed field for small-scale fluctuation dynamo. 
The asymptotic exponential growth rate and the small- 
scale magnetic correlations are exactly the same as for 
any other random seed field whose correlation function is 



non-orthogonal to the leading dynamo mode. We there- 
fore do not agree with authors [33-35 who distinguish 
between "fluctuation dynamo" and "magnetic induction" 
as two fundamentally different mechanisms. 

For example, Schiissler and Vogler [35] have proposed 
to explain small-scale magnetic fields in the quiet solar 
photosphere as a consequence of near-surface dynamo ac- 
tion. They mention magnetic induction as a possible al- 
ternative explanation, but dismiss it with the remark: 
" 'Shredding' of pre-existing magnetic flux (remnants 
of bipolar magnetic regions) cannot explain the large 
amount of observed horizontal flux since the turbulent 
cascade does not lead to an accumulation of energy (and 
generation of a spectral maximum) at small scales. On 
the other hand, such a behavior is typical for turbulent 
dynamo action." We find in the KK model that, quite 
to the contrary, induction of a mean field produces the 
same small-scale fluctuations and energy spectra as does 
the fluctuation dynamo. We may indeed say that — at 
least in the kinematic regime of weak fields — "magnetic 
induction" is nothing but "fluctuation dynamo" with a 
large-scale, deterministic (mean) seed field. 

Although Schekochihin et al. [33] do make a distinction 
between two separate mechanisms, their final conclusion 
is not so different from ours. In their Fig. 8(a) they find, 

(c) 

in a sub-threshold regime with Re m < Re m \ that satu- 
rated energy spectra for induction from a uniform field 
and for decaying spectra of the "failed" dynamo state ex- 
actly coincide, when normalized to the same total energy. 
They conclude that "the same mechanism is responsible 
for setting the shape of the spectrum of the magnetic 
fluctuations induced by a mean field and of the decaying 
or growing such fluctuations in the absence of a mean 
field," just as we see in the KK model. 

We have also studied the KK model with velocity 
roughness exponent £ < 1, where the Pr m = fluctua- 
tion dynamo does not exist, and with a uniform magnetic 
seed field in order to get some insight into the physics of 
the induction mechanism in a failed-dyamo regime. This 
is very roughly the same situation that was studied in sev- 
eral liquid-sodium experiments with Pr m < 10 -5 [27U3T] 
and in related numerical simulations [32 , 33 at somewhat 
larger Prandtl numbers. Some of those studies have re- 
ported observing a k~ 5 / 3 spectrum of magnetic fluctua- 
tions in the velocity inertial-range [29J [31j [32] , while oth- 
ers have reported a k~ x spectrum [28[|33j. The very least 
we can conclude from our analysis is that the explana- 
tions offered for those spectra based on inertial-range tur- 
bulence physics do not successfully explain the results in 
the KK model, even where those explanations should ap- 
parently apply. An exception is the Golitsyn-Moffatt ar- 
gument for the spectrum at wavenumbers k > k^, which 
succeeds (with slight modification) in the KK model. At 
intermediate wavenumbers /cz,(t) < k < k^, we find in 
the KK model a range with input of magnetic energy by 
induction of the mean field balanced by nonlinear stretch- 
ing and cascade to the resistive scale. This is the same 
physics invoked in the theory of Ruzmaikin and Shukurov 
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[44] , who predicted a spectrum (Bfyk -1 on dimensional 
grounds. However, this prediction is only verified in the 
KK model for £ = 0, where nonlinear cascade can be 
accurately modelled as an eddy-diffusivity. For the KK 
model with < £ < 1 this prediction is modified by a 
large anomalous exponent (with |£i| = 2 for £ = 1) 
which must be calculated by a non-perturbative argu- 
ment. At low wavenumbers k < /cl(£) in the KK model 
we find a range with continuous growth of magnetic en- 
ergy spectrum supplied by induction of the uniform field. 
Unlike the argument in [29] [3TJ [32] , there is no balance 
with large-scale advection, although such effects exist in 
the model. Indeed, it is very unclear to us how large-scale 
advection, which conserves magnetic energy, can by itself 
balance the input from induction. 

We shall not attempt here to explain in detail the in- 
duced spectra observed in the experiments and simula- 
tions, except to note that they must involve large-scale 
physics outside inertial-range scales. This was, in fact, 
the point of view of Bourgoin et al. [28] who explained 
their k~ x spectrum by global fluctuations of the flow pat- 
tern. Large-scale effects are certainly necessary to obtain 
saturation of the energy. If the length-scale L(t) in (74) 
had a finite limit for large times, e.g. for turbulence 
confined to a box of size L^, then the system should 
reach a statistical steady- state, just as in the randomly- 
forced case [10]. Otherwise, magnetic energy would con- 
tinue to grow without bound within a kinematic descrip- 
tion. In such a steady-state, the only extended spectral 
ranges would be the nonlinear range ( 74 ) and the Golit- 
syn range ([73|. The latter is well-confirmed in experi- 
ments [27] [28J [3T] and simulations [32[ [33]. If the k~ x 
or & -5//3 spectra observed at larger scales can be under- 
stood at all in terms of inertial-range turbulence, then it 
seems that an analogue of the nonlinear range ( 74 ) is the 



most likely explanation. If so, then the determination of 
the precise exponent is a hard theoretical problem, for 
there may be large anomalous scaling corrections to the 
dimensional Ruzmaikin-Shukurov k~ x spectrum. 

The precise spectral exponent of induced magnetic 
fluctuations in the liquid metal experiments is open to 
debate. Over the limited scaling ranges available it is 
quite difficult to distinguish between —5/3 and —1 expo- 
nents and the empirical spectra can be fit about equally 
well with both power-laws, or others as well. It is not 
even clear that the concept of "scaling exponent" is en- 
tirely well-defined without the possibility to extend the 
putative scaling ranges. Our three spectral ranges (73), 
(74),(75) can be made arbitrarily long by adjusting pa- 
rameters, but, in the experiments and simulations, only 
the Golitsyn spectral range can be lengthened by low- 
ering Pr m . The low-wavenumber spectral ranges can be 
made longer only by increasing Re mi but, beyond a criti- 

(c) 

cal value Rem , the dynamo effect sets in and the physical 
phenomena change. 



Appendix A: Boldyrev-Cattaneo-Rosner formalism 

We consider in this appendix kinematical relations be- 
tween the two-point correlations of magnetic fields and 
vector potentials, first for spatially homogeneous ensem- 
bles and then for both homogeneous and isotropic en- 
sembles (but not reflection-symmetric). The results are 
closely related to those in classical references on homoge- 
neous turbulence, such as [46^ and [47 , but we employ a 
mathematical reformulation of [14 (hereafter, BCR) that 
was developed for helical turbulence. We discuss here, in 
particular, the consequences of gauge- invar iance. 

The two-point correlation of the magnetic field in 
a homogeneous ensemble is defined as C u '(r,t) = 
(B l (r,t)B^(0,t)). This satisfies the solenoidality condi- 
tions diC l i = djC l i = 0. The two-point correlation 
of the magnetic vector-potential is likewise Qij (r, t) = 
(Ai(r, i)Aj(Q, t)). A pure gauge field A\ = di\ has corre- 
lation Q\. = -didjA with A(r) = (A(r)A(O)). Solenoidal 
correlations C and pure-gauge correlations Q 9 form mutu- 
ally orthogonal subspaces in the Hilbert space with inner 
product 



(C 



tf'rC«(r)ey(r). 



(Al) 



This observation will prove important in what follows. 
Notice also that the magnetic and vector-potential fields 
are related by curls, so that 

C ij = _ V ioMg M 

with the non-positive, self-adjoint differential operator 
v ijM = e ^P e o^d p d q . 

It is furthermore useful to introduce the joint correlations 
of magnetic fields and vector-potentials 

* 4 fc (r, t) = {B\v, t)A k (0, t)>, * fc 4 (r, t) = (A k (r, t)B*(0, t)> 

which are related by fc (r, t) = ^/ k % (— r, t). These are 
obtained from the vector-potential correlation by 



with 



and likewise 



n i,k e ipk g 



'pi 



The symbol * here denotes adjoint with respect to 
(f,g) = fd 3 rf(r)g(r). Note that V^ ke = -K^iTV* 1 )*. 

We are mainly concerned with statistics both homoge- 
neous and isotropic, but not necessarily invariant under 
space reflections. In that case, the magnetic correlation 
becomes 



C ij (r,t) = C L r l P + C N (5 ij - r l r J ) + C H e ljk r k , 
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with coefficients depending only on r, t, and the 
solenoidal condition becomes 

C N (r,t)=C L (r,r) + ±rC' L (r,t). 

Note that there is no condition on C#, as the third term 
is always divergence- free. One can similarly write the 
vector-potential correlation as 

Gij(r,t) = GLViVj + G N (Sij - fiTj) + G H e ijk r k . 

For a pure-gauge field 

G 9 L (r, t) = A"(r, i), G%(r, t) = *A'(r, t), G a H {r, t) = 0. 

This is equivalent to the conditions 

GlM) = G N (r,t) + rG'^r,*), G H (r,*) = 0. 

We consider also the mixed magnetic and vector po- 
tential correlators: 

¥ ^r, *) = tt^f + ^ N (S^ - r'rj) + ^6%. fc f fe , 

This correlation must be divergence-free in both indices i 
and j. This is trivial for i and follows in j for the first two 
terms by their symmetry in z, j. The third antisymmet- 
ric term is always solenoidal in both indices. Another 
interesting fact is that this mixed correlation is gauge- 
invariant in the isotropic sector. To see this, note that 

<B*(r,t)A((M)> = P(r,ty 

by isotropy. The solenoidal condition implies that 
P(r,i) = a(t)/r 2 and regularity for r — >> implies that 
a(t) = 0. By differentiation it then follows that also 

(S < (r,t)5 fe A(0,t)) = 0, 

which implies gauge-invariance of ^ l k . The same argu- 
ment implies also that every vector-potential correlation 
in the isotropic sector can be written as a sum Q = Q Sj rQ 9 
of a solenoidal correlation Q s and pure-gauge correlation 
Q 9 . This follows by writing the vector-potential field as 
A = A s + VA, where A s is the Coloumb-gauge potential 
satisfying V*A S = 0, and Q s is thereby identified as the 
correlator of the Coloumb-gauge vector potential field. 

The relation & l £ = lZ % ' k Q k g becomes in the isotropic 
sector 



* L = 1—, * N = — + G' H , ^ H = 



Gl — Gn 

r r " ' r 



-G' r 



N ' 



The solenoidality condition S&n = + \ v ^'l ls satis- 
fied automatically, consistent with our earlier conclusion. 
The relation = lZ % ' k ^ k likewise becomes 



C L , = 2-^, C N = ^ + *' H , C H = * L * N -*' N 



r 



by using 

(r, t) = ^ L r k r j + ^ at(V ^) + ^He k j£ f e , 



with the same coefficient scalar functions as Note 
also that C H = -2^' L - by employing the 

solenoidality relation between S&n and ^z,- 

The main idea of BCR was to employ a slightly differ- 
ent basis, 



C i \v,t) = C 1 



e ijk r k 



— — h C 2 h C 3 — p- 

V^r r y/2r 



and likewise for Qij, with normalizations chosen so that 
the Hilbert space inner product in the isotropic sector 
becomes 

(C, g)=Aix dr [CiGi + C 2 G 2 + C 3 G 3 ]. 
Jo 

Trivially d = V2rC N , C 2 = rC L , C 3 = V2rC H - The 
main advantage of this normalization is the simplifica- 
tion of various differential operators. A straightforward 
calculation (see details below) gives C = — VQ with 










acting on the 3- vector function Q = (Gi,G 2 ,G3) T . The 
operator V is explicitly self- adjoint in this representation. 
Similar simplifications occur in the operators 1Z and TZ* . 
Following BCR we introduce 



W 



with 

W 2 = ^3 = Vlr^H, W s = ^1 - V2^ 2 = -^rV L . 

This representation exploits the fact that \I/ has implic- 
itly only two degrees of freedom, due to the solenoidality 
condition on S&l and ^ n- A bit of calculation yields 



with 



iz 




c = nw, w = n*g 




The range of the operator IZ is easily checked to consist 
of the solenoidal correlations, while its adjoint 7Z* anni- 
hilates pure-gauge fields, associated to gauge-invariance 
of W. Combining the above results gives the previous 
formula for V = —1Z1Z*. We see that V is a negative- 
definite, self- adjoint operator in the isotropic sector. 
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Appendix B: Numerical Methods 

This appendix discusses the numerical solution of the 
equation (55) in the text: 



d r (r 2 bd r A) + Ar 2 A = -r 2 F. 

The output of MATSLISE provides an input to this equa- 
tion, via the source term F = 2$ + TG. We shall there- 
fore employ the same grid of points r^, i = 1, that 
was generated by the MATSLISE algorithm for solving the 
Sturm-Liouville problem (44). It is useful to observe that 



that grid is approximately evenly spaced on a logarithmic 
scale. Thus, we introduce into ( [55] ) the variable a = lnr 
and multiply by r, giving 



d CT (p(a)d CT A) + q(a)A = -f(a) 
with definitions (different than in ([44])) 



(Bl) 



p = e°b{e°), q = \e 3 °, f 
We then discretize ( |B1[ ) as 
'A i+i - A 



3a 



1 

8 a i 



S(Ti 



Aj ~ Ai-i 



with 8(Ji = <Ji+i — (Ji and G{ = lnr^, i = 1,2,...,L . An 
extra point is added with ol+i = oo and A^ + i = to 
define the components with i or j = L. This leads to the 
L- vector equation 

MA = F (B2) 
with negative-definite, self- adjoint, tridiagonal matrix 

ilfe = — + — + 



&7i 



q{(Ti)0(Ti 1 & 



and vector F 7 - 



Scii da i-\ 
f{(7i)8(7i. The z = L equation is 



integral J Q r dp a(p)W2(p) / p using the composite trape- 
zoidal rule (cumtrapz in MATLAB). The first (L — 1) linear 
equations from ( |B2[ ) are then solved using mldivide in 
MATLAB, giving A*, i = 1, 2, L - 1. 

Getting corrections 8G X = y/2A'(r), 8G 2 (r) = rA"(r) 
requires us to approximate also derivatives of A(r). We 
employ the formula A'{r) = \d a A discretized using cen- 
tral differences as 



A ■ = exp 



A 



■t+i 



A, 



^z+l 



(Ji 



At the endpoints i = 1,L corresponding forward and 
backward-difference approximations are employed. To 
calculate the second-derivative A^, the same formulas are 
employed with A^ replaced by A^ . The difference approxi- 
mations yields some spurious oscillations for r < 2, which 
are removed by filtering. 

The parameters Aq,Ai, and A s in the small-r asymp- 
totics (56),(57) for A(r) were determined as follows. We 
found directly that A = —1535.5 by taking the limit of 
small r in the numerical solution. We then calculated 
Ai in two different ways. First, we solved eq.(58), or 
AAo + 12Ai = — Fq. To get a sufficiently accurate value 
of Fq = 2<3>o + 24y/2A/X we had to calculate an accurate 
value of the integral 



dr- 



aW 2 



-208.1745 



To get this value, we extended the numerical solution for 
W2 from MATSLISE using the large-r asymptotic formula 
(47) from r = 560 to r = 10 4 , and then used a compos- 



ite trapezoid rule to approximate the integral. This gave 
F = 439.043 and a resulting value Ai = 11.9299. On 
the other hand, we could also determine Ai by a least- 
square linear fit to A(r) — Ao in a log- log plot. Over the 
range from r = 2.7 x 10 -6 to r = 0.12 the best fit was 
A(r) - A = Air m with m = 1.9999 and Ai = 11.9517, 
agreeing quite well with the previous, independent deter- 
mination. Finally, we calculated A s = —6.0232 from the 



already determined numerical constants and eq.(59). It 



was verified that the full small-r asymptotic expression 
A = A reg + A sing with these constants gave an accurate 
fit to the numerical result for r < 0.5. 



q(cr L )A L = f(a L ) 

with direct solution A^ = F(cr L )/A. The remaining L — l 
equations can be solved with the above value as bound- 
ary condition on A^ or, more simply, with the boundary 
condition A^ = (since F(cr) —> as a — >• 00). We have 
verified that both choices of boundary conditions lead to 
very quantitatively similar results; we thus present only 
those with the condition A^ = 0. To calculate the quan- 
tity <I> that appears in the source-term F we evaluate the 
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